Poor man's derivation of the Bethe-Ansatz equations for the Dicke model 
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We present an elementary derivation of the exact solution (Bethe-Ansatz equations) of the Dicke 
model, using only commutation relations and an informed Ansatz for the structure of its eigenstates. 



In 1954, Dicke showed that a model describing a set of 
two-level systems coupled to a quantised electromagnetic 
mode leads to a supperradicnt effect^. Generalisations to 
multicomponent systems naturally appear in various ex- 
perimentally relevant contexts 2 ^. Other generalisations 
involve a spatially extended photonic field^, or itinerant 
two level systems^, motivated by an experiment on cold 
atoms in a 2D lattice coupled to an optical resonator—. 

The Bethe-Ansatz solution for the Dicke model with in- 
homogeneous excitation energies was originally obtained 
by Gaudini as a side result of solving the central spin 
problem. Using a variational method that results in 
complex algebraic computations, he showed that the so- 
lution of the central spin problem is equivalent to the 
Bethe-Ansatz solution of the BCS problem derived by 
Richardson^. By expanding the Bethe-Ansatz equations 
for the central spin model in the limit of large central 
spin, Gaudin obtained corresponding equations for the 
Dicke modeli. Though this procedure solves the original 
problem, the derivation is computationally complex, and 
thus not easily extended to other, related models. 

The purpose of this paper is to provide an elementary 
derivation that starts from the original Dicke model, in 
the hope that our simplified treatment might pave the 
way toward finding similar solutions to generalized Dicke 
models. We follow a method suggested by Richardson-, 
for the BCS model and presented in Ref. [TolfTlI This 
methods exploits the observation that the structure of 
the exact eigenstates of the Dicke model is similar to 
that of an auxiliary model, involving only bosons. The 
only difference is that the eigenvalue equations that de- 
termine the quasi-energies characterizing these states be- 
come more complicated for the Dicke model: they turn 
into Gaudin's Bethe-Ansatz equations, which we derive 
here using only commutation relations. 

The inhomogeneous Dicke model describes a set of 
non-identical two-level systems with excitation energies 
Cj and a single photon mode with frequency u, coupled 
with interaction strength g: 

N 

H = cub^b + £ e 3 {SI + \) + g J2 (Sp + Sjtf) . (1) 

3=1 3 

The spin-i operators satisfy (S^) 2 = and 

& , S+] = -2S*6 t] , [5?, Sf] = ±SfSii , (2) 
while the boson operators satisfy \b,b'] = 1. 



Let | Vac) be the "vacuum" state containing no boson 
excitations and all spins down, i.e. b\ Vac) = S~ | Vac) = 

0. H commutes with the operator b^b + y^j—i Sji which 
counts the number of excitations relative to |Vac). Thus, 
P-eigenstates can be constructed by acting on |Vac) with 
(products of) linear combinations of Sj~ and b^ operators, 
of the general (unnormalized) form 

N 
3=1 

where the coefficients A V j arc to be determined. For an 
cigenstate with n excitations relative to |Vac) we thus 
make the Ansatz (followin g 10 ' 11 ), 

|tt„> =/T|Vac) , (4) 

where we use the shorthand notation (for n' < n) 

n 

p r :> = n B t (5) 

v— n' 

for a product of B^'s (for n' > n, we set P™/ = 1). For 
later use, note that such products satisfy the composition 
rule P%P? +1 = P% for n' < v < n. 

Wc require that H\^ n ) = £ n \^ n ). Commuting H past 
P™ to the right and using 7J|Vac) = 0, wc obtain 

(£„P 1 "-[i?,P 1 l ])|Vac) = 0. (6) 

Using the general operator identity 

n 

[x,p:,]=y j p^ 1 \x,bI]p: +1 , (?) 

u—n' 

Eq. ^ can be written as 

(^P? ~ J2 Pr\H, Pt]p» +1 j | vac) = . (8) 
The requisite commutator is given by 

N 

[H, Bt] = J2( A ^3 + 9)S+ + (w - 2gX v )b\ (9) 

3=1 

where X v = S^Li A V jS*. By making the choice 

A» 3 = ~~p~~~~~7~ , (10) 
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where the parameters E v will be called quasi-energies, 
Eq. §§§ can be brought into the simplified form 

[H, Bt] = E v Bl + {u-E v - 2gX v )tf . (11) 

Inserting this into Eq. (jSJ and identifying the eigenergy 
with the sum on quasi-energies, £ n — Y] , E v , yields 

n 

Pr\" -E v - 2gX v )P? +1 tf | Vac) = . (12) 

To make sense of this condition consider, for a moment, 
an auxiliary, purely bosonic model, obtained from the 
Dicke Hamiltonian ((T|) by replacing Sf, Sj~ and (SJ + |) 

by frt, bj and bjbj, respectively, with [6^, ftt] = 5%. Re- 
peating the above analysis yields only one change: since 
[6j , bt] gives 1 instead of [S~ , Sj~] giving —2S?, the op- 
erator X v in Eq. (|9|) is replaced by the c-number x v = 
-ISiLiAtf. Thus Eq. (Hl|) can be satisfied by requir- 
ing that lo — E v — 2gx v = for all v. Via Eq. (|TU|) this 

implies ui — E v + J2j=i 9 1 l{P*v — e j) = 0; which deter- 
mines the E v . This equation can also be obtained by 
making the Ansatz H = Y)„ E V B\,B V and demanding 
that [H,Bl] = E V B\. For this auxiliary model the B\, 
thus describe independent single-particle excitations, and 
the quasi-energies E v are their eigenergies. 

Let us now return to the Dicke model, where X v is an 
operator, so that we have to work a little (but not much!) 
harder to satisfy Eq. (fT2"j) . To this end, commute X v past 
P"_i_i to the right and use -X"„|Vac) = cc„|Vac), to obtain 



E Pr lp ?+i(" - E »- 2gx„p\ Vac) (13a) 

n 

= 2gY J Pr 1 \X u ,Pl\ 1 p\\ & c) . (13b) 

To simplify the second line, use Eq. ((7]) and the relation 

Bl-Bl 



to write E u =iPi~ [^.C+i] as 



(14) 



" ™ Bt _ Rt 

-gY.p'r 1 E /^Vi (is 
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E^r x c +1 E 
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E v — En 



(15b) 



Eq. (|15b[) follows by relabelling ^ o /i in the Bt term 
of Eq. (fT5aj) . Inserting Eq. (|15b|) into Eq. (|13b|) . we 
note that Eq. (fT3"|) is satisfied provided that the n quasi- 
energies E v obey the following n coupled equations: 



E v 



N 
3 = 1 



E ■ 

fi—l,fiy^iy 



2ff 2 



(16) 



which follows from A V jA^ = —g(A U j — A^j)/(E V — Ey), 



These are the celebrated Bethe- Ansatz equations for the 
Dicke model, first obtained by Gaudini. The fact that 
the right-hand side couples the equations for different E v 
together presents the additional complication arising for 
the Dicke model in comparison to the above-mentioned 
auxiliary boson model. It implies that the B\, do not 
describe independent single-particle excitations, since the 
value of any E v depends on that of all others. 

Generally Eqs. (fH)]) have to be solved numerically. 
For sufficiently small n, however, the original model 
(P) can be diagonalised directly by solving the eigen- 
value problem in the basis of uncoupled bosonic and spin 
cigenstatesi^ instead of the basis (0}. 

It is straightforward to expand the normalization fac- 
tors of Gaudin eigenstatesi and verify that | (^ n \^n) | 2 = 
det M, where M is an n x n matrix with elements 

M w = i+Ef =1 ^.-2E: =lli 

and we used the shorthand A vfl = gj (E v — E^). 
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